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Abstract 

We investigate the associated production of a scalar Higgs boson (h° or H°) with Z° boson in 
the minimal supersymmetric extension of the standard model (MSSM) at the CERN Large Hadron 
Collider (LHC), including the contributions from bb annihilation at the tree level and gluon fusion 
via quark and squark loops. We quantitatively analyze the total cross sections in the mSUGRA 
scenario. For the production of hP associated with Z°, we find that in most of the parameter regions, 
the contributions from initial bb and gg are at a level of one percent of the total cross section and 
therefore almost insignificant. For the production of H° associated with Z°, the contributions from 
bb channel can be much larger than those from light quark initial states. Especially for large tan f3, 
the increment can reach about one order of magnitude. Thus, when considering the associated 
production of H° and Z° at the LHC, the contributions from bb annihilation should be taken into 
account seriously. 
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I. INTRODUCTION 



The search for Higgs bosons is one of the main goals of the CERN Large Hadron Collider 
(LHC), with y/s = 14TeV and a luminosity of 100 fb _1 per year 4]. In the standard 
model (SM), the Higgs boson mass is basically a free parameter with an upper bound of 
m H < 600 — 800GeV 2]. However, present data from precision measurements of electroweak 
quantities indicate the existence of a light Higgs boson (m# < 204 GeV at 95% C.L.) and 
direct searches rule out the case m# < 114GeV [3|. In addition, in various extensions of the 
SM, for example, in the two-Higgs-doublet models (THDM) [4], particularly in the minimal 
super symmetric standard model (MSSM) t 5], there are five physical Higgs particles: two 
neutral CP-even bosons h° and H°, one neutral CP-odd boson A , and two charged bosons 
H ± . The Higgs boson h° should be lightest, with a mass m^o < 135GeV when including the 
radiative corrections [6j . It has been shown [7( that the h° boson cannot escape detection at 
the LHC and that in large areas of the parameter space, more than one Higgs particle can 
be found. 

At the LHC, the neutral Higgs bosons can be produced through following mechanisms: 

gluon fusion gg — > </> B fl [J Q, weak boson fusion qq - qqV*V* - qqh°/q q H° Q, 

I I 

associated production with weak bosons [13], associated production with a heavy quark- 
antiquark pair gg,qq —>■ ti<p/bb(p Q] and pair production Q]. In this paper we focus our 
attentions on the production of h°/H° in association with Z° boson in the MSSM. This 
is one of the main discovery channels of neutral Higgs bosons at the Tevatron, with the 
Higgs decays to bb and Z° decays leptonicly. The main backgrounds to this signal are Zbb 
production, ZZ production and top quark production [lj|. At the LHC, it can be used to 
detect the "invisible" Higgs which decays into the LSPs. The detailed analysis of signals and 
backgrounds can be found in Ref. IT], where they found the ratio of signal to background 
can as large as 12. In the SM, the Drell-Yan production process (light qq annihilation) of 
Higgs associated with Z° boson [l3| . as well as the gluon fusion process jla], has been studied 
previously. In the MSSM, the Drell-Yan contributions are related to the ones in the SM 
by an overall coefficient, but the contributions of gluon fusion are generally different in the 
two models due to the different couplings and additional contributions from new channels. 
Meanwhile, there are potentially important contributions to the production of neutral Higgs 
bosons from bb annihilation at the tree-level. In the MSSM, the Yukawa couplings b-b-cf) 
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can enhance Higgs boson production cross section via bb annihilation significantly for large 
tan f3. Therefore, besides the production channel via Drell-Yan process, in order to obtain 
the complete leading order (LO) total cross sections for the associated production of scalar 
Higgs bosons with Z° boson, the partonic subprocesses bb — > Zh°/ZH° should also be taken 
into account. Moreover, the loop induced subprocess gg — > Zh°/ZH° become important 
due to the potential enhancement, and should be considered, too. We will investigate all 
these contributions in our work. 

Our paper is organized as follows. In section |Hl we list analytical results for the tree 
level cross sections of pp — > bb — > Zh° /ZH° and loop induced cross sections of pp — > gg — > 
Zh°/ZH° in the MSSM. In section UTT1 we present quantitative predictions for the inclusive 
cross section of pp — > Zh°/ZH° + X at the LHC adopting the MSSM parameters constrained 
within the minimal supergravity (mSUGRA) scenario and discuss the implications of our 
results. The relevant MSSM couplings and form factors are given in the Appendix [X] and 
IB] respectively. 

II. CALCULATION 

I — I 

The relevant Feynman diagrams are created by FeynArts |l9| version 3.2 automatically 
and are shown in Fig. We carry out the calculation in the 't Hooft-Feynman gauge 
and use dimensional reduction for regularization of the ultraviolet divergences in the loop 
diagrams. In the following expressions, G^ R and G^ k are the couplings, which are given 
explicitly in Appendix [XJ H stands for the scalar Higgs bosons, h° or H°; S and / are 
propagating scalar and fermion particles, respectively. 

For the partonic subprocesses 

b(h) + b(k 2 ) -> Z(k 3 ,e 3 ) + H{h), 
g(ki,e 1 )+g(k 2 ,e 2 ) -> Z(k 3 , e 3 ) + #(fc 4 ), 

we define the Mandelstam variables as 

s = (h + k 2 ) 2 , t = (k l -k 3 ) 2 , u = (k 1 -k 4 ) 2 . (l) 

The tree- level amplitude M. bb for the subprocess bb —> ZH consists of the three diagrams 
(a) — (c) in Fig. [TJ The previous results [13j] for vector boson bremsstrahlung include only the 
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first diagram, where the bb pair replaced by light quark- ant iquark pairs. We also recompute 
the contributions of light quarks to the tree-level cross sections and compare them with ones 
of bb -> ZH. 

Using the notations defined above, the amplitude Ai bb can be expressed as 



M bb = M bb + M bb + M bb ( 2 ) 



with 



-iG HZZ 

M bb = v(k 2 ) MGTPl + G bbZ P R )u(h) - _ 2 , (3) 

z 



s — m 



lCr HSZ 

M bb = J2 vik^G^Mhh TT . F (h + 2fc 4 ) ■ e 3 , (4) 

S=ao,G° S ~ m S + lrn S*- S 

M c = v{k 2 )G j M&l P L + G R P R )u{ki) 

t — mf 

+ v(k 2 ) MG^Pl + G lbZ P R ) l{h Z A t mb) g^(fei)- (5) 

u — mi 

Here Pl,r = (1 =F 7 5 )/2, mgo = m^, r G o = 0, and T^o is the decay width of A . 

The gluon fusion subprocess is forbidden at tree-level. At one-loop level, in general, the 
cross section will receive contributions from both quark loops and squark loops, as shown 
in Fig. |2J El Note that each diagram actually represents a couple of diagrams with opposite 
directions of charge flow. However, we find that the contributions from each pair of squark 
loop diagrams to this process cancel each other due to the opposite signs of momenta. So 
the gluon fusion cross section arises only from the quark loop diagrams, i.e., 

M" = M 99 + M 9 b 9 + ■■■ + M 99 , (6) 

where the subscripts a — e refer to the corresponding diagrams in Fig. |21 

In general, the amplitudes M." can be written as a linear combination of the invari- 
ants formed by three independent external momenta ki,k 2 ,ks, polarization vectors of three 
external gauge bosons £i,e 2 ,Ez, metric tensor g^ v and Levi-Civita tensor e^ vpa , in which 
the terms without Levi-Civita tensor vanish. Thus, with taking into account the relations 
E\ ■ ki = e 2 ■ k 2 = £3 ■ k 3 = 0, 24 terms remain in the amplitudes. It is easy to prove the 
following identity 

gHU e paa/3 _ gHp e uaa/3 _ ^pa^vpafi _|_ gpf3 ^upoa ^ 
_|_ g^°" e MP"/3 _ gpo e pvafi _|_ gaa e pupl3 _ ^j^^vpa _ g 



Using it, one can immediately see that not all of the 24 terms are linear independent. 
Actually, there are only 14 independent ones, and finally the amplitudes can be written as 

M" = A ie £l£2£3/Cl + A 2 e £l£2£sk2 + A 3 e £l£2£3ks + A A e £l£2klk2 e 3 ■ h + A 5 e £l£2klk2 e 3 ■ k 2 
+A 6 e £l£3klk3 e 2 ■ h + A 7 e £l£sklks e 2 ■ h + A 8 e £2£3k2k3 e l ■ k 2 + A 9 e £2£3k2ka e x ■ k 3 
+A 10 e £lklk2k3 e 2 ■ e 3 + A^e^^Sx ■ e 3 + A l2 e £sklk2ks e x ■ e 2 (8) 
+A 13 e £l£2k2k3 e 3 ■ h + A 1A e £l£2klk3 e 3 ■ k 2 . 

The explicit expressions of the form factors Ai(i = 1, 14) are shown in Appendix [Bj 
The differential cross sections of the subprocesses are given by 

da 1 



■Ei^r. ( 9 ) 



dt 16tts(s — Am 2 ) ■ 

where M. = A4 bb , m = rut, for bb channel and Ai = A4", m = for gluon fusion. 

Using the standard factorization procedure, the total cross section can be obtained as a 
convolution of partonic cross sections with corresponding parton distribution functions, 

a(pp -> ZH) = i , r / dx! / dx 2 fa/pMfp/pM + [a «-> (3) cr(a(3 -> ZH), 



(10) 

where r = (mz + m//) 2 /s, y/s is the center-of-mass energy of the LHC, a and j3 denote the 
initial partons. 

III. NUMERICAL RESULTS AND CONCLUSIONS 

In the following we present some numerical results. In our numerical calculations, the SM 
input parameters were taken to be a(mz) = 1/128.8, m w = 80.423GeV, m z = 91.188GeV 
and rru = 174.3GeV [3]. We use the one-loop evolution of the strong coupling constant a s (Q) 
2j], |21( with a s (mz) = 0.1172. The leading order CTEQ6 parton distribution functions 
221 are used here and the factorization scale is taken to be the invariant mass of the two 



final particles fj, = m ZH = yipz + Pn) 2 - Moreover, in order to improve the perturbative 
calculations, we take the running mass of bottom quark m^Q) evaluated by the one- loop 
formula 

m b (Q) = U 6 (Q,m t )U 5 (m t ,m b )m b (m b ), (11) 
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with mb{rrib) = 4.25GeV. The evolution factor Uf is 

a 

as{Qi) 



uAQiM = (^M)' (12) 



with 

d > = whf- < 13 > 

The relevant MSSM parameters, the Higgs masses and mixing angles a, are determined 



in the mSUGRA scenario as implemented in program package ISA JET 7.69 |23j]. The 
GUT parameters m , A and sgn(//) were taken to be m = 200 GeV, A = —100 GeV, 
fi > 0. mi/2, tan/5 were varied to obtain various Higgs masses. Actually, in mSUGRA 
scenario, the mass of the light scalar Higgs boson h° can only reach about 125 GeV, which 
is independent of the choice of GUT parameters. It should be noted that in some parameter 
region, m^o > mz + ttt^o, and the momentum of A can approach its mass shell, which will 
lead to a singularity arising from the A propagator. This can be avoided by introduing the 
non-zero decay width T^o, which was also calculated by ISA JET. 

Fig. Wi a ) and Efa) show the total cross sections for the process pp — ► Zh° + X versus 
the light scalar Higgs boson mass m^o for tan/3 = 4, 15 and 40. Fig. 0^6) andE^fr) show the 
dependence of the cross sections on tan (3 for m^o = 105 and 115 GeV. From Fig. H] we can 
see that the bb contributions are approximately one order of magnitude smaller than the 
ones of Drell-Yan process ,13|, and only increase the total cross sections by about several 
percents, which are smaller than the QCD corrections to the Drell-Yan process. Fig. |3] 
shows the contributions of gluon fusion to the cross sections. Comparing with ones of the 
complete qq annihilation (i.e. Drell-Yan + bb), one can see that for a light h (for example, 
rriho = 105 GeV), the contributions of gluon fusion are about half of the ones of the qq 
annihilation for low tan/3(< 5), but with the increasing of tan/3 from 4 to 10, the former 
decrease significantly, and become one order of magnitude smaller than the latter in general. 
Therefore, in most of the parameter regions, the contributions of both bb and gg initial states 
are almost insignificant and can be neglected. 

Fig. E(a) and Ufa) show the total cross sections for the process pp — ► ZH° + X as a 
function of the heavy scalar Higgs boson mass m#o for the three representative values of 
tan f3. Fig. WLp) and Ulp) show the dependence of the cross sections on tan f3 for m#o = 200 
and 400 GeV. We found that the bb contributions to the total cross sections for the production 
of heavier neutral Higgs boson can be much larger than the contributions of Drell-Yan process 



for large tan/3, as shown in Fig. El For example, when tan/3 = 40 and m#o = 200 GeV, the 
bb cross section is about 50 fb, while one of Drell-Yan process is less than 0.5 fb. Therefore, 
the bb channel provides a dramatic enhancement to the total cross sections, and should 
definitely be considered. On the other hand, the contributions of gluon fusion are very small 
(< 1 fb) and can be neglected in most shown in Fig. [7] 

For the comparison of the production rates of H$m, h° and H°, they are displayed as 
the functions of their masses in Fig. IS] where all the contributions (i.e. Drell-Yan + bb + 
gluon fusion) have been included for tan/3 = 4 and 40, respectively. For the SM Higgs, as 
shown in Ref. Q|, the gluon fusion contributions are negligible for Higgs masses allowed 
by experiments, so we didn't include them here. From Fig. |SJ one can find that no matter 
what the value of tan /3 is, the production rates of H° are the smallest, the ones of h° are the 
largest, and the ones of H$m are medium. This feature indicates that the predictions for 
the associated production of the Higgs bosons and Z° boson in the SM and the MSSM are 
different quantitatively and distinguishable, which is in agreement with the results shown in 
Ref. [2l| . where only Drell-Yan contributions are included. 

In conclusion, we have calculated the scalar Higgs bosons (h°, H°) production in associ- 
ation with a Z° boson through both bb channel and gluon fusion in the mSUGRA model 
at the LHC Our results show that in most of the parameter regions, the contributions to 
the total cross section for associated production of the light scalar Higgs boson h° and Z° 
boson mainly come from light quark annihilation and the contributions from initial bb and 
gg are small and negligible. And the total cross section for associated production of heavy 
scalar Higgs boson H° and Z° boson, with including the contributions of bb channel, can be 
increased greatly. Especially for large tan/3, such increment can reach about one order of 
magnitude. Thus, the contributions of bb channel should definitely be taken into account, 



171 ] . which is based on the previous 



and the ratio of signal to background obtained in Ref. 
calculations, can be enhanced at large tan/3. The contributions from gluon fusion are still 
small and unimportant. 
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APPENDIX A: COUPLINGS 



Here we list the relevant couplings in the amplitudes, 
and r, s stand for color indices. 

ieS ij 



QV?V? Z 



Qd % (PZ 



- 3 + 4 4)t^— > G T 3Z 



i,j stand for generation indices 



2ies w 5 lj 



Wj 



u^h 



-(-3 + 2s 

ieM w S0. 

ic a eS^m u i 



Qcwsw.' 



6c w s w 



2M w s@s w 
„j,r ies a 5 tJ m ul 



G 

G 

G 

G 

G 
G 
G 



u^G 



2MwSpSw 

eb % im u i 

2M\yS\ytf3 



2M w s w ' 
h°A°z _ c p-a e 

2cyySyy 
H°A°Z _ eSf3- a 



r>d % dPZ 



qH u ZZ 

Qd^iH 
d^d^A 



G 



G 



d^^G 



3cw 

3c w 
ieM w cp_ a 

2 ' 
C w S\v 

2cpM w s w ' 

2cpM w sw 
etp5 l3 m d i 

2M w sw ' 



G 



h°G°Z 



H°G°Z 



G 

Qd^dlg" 



2M W S W 

eS/3-g 

2cyySw 

cp- a e 

2cy/Sw 

-ig s S ij T r a s . 



with 



sw = sin ^vy, cw = cos9w, s a = sin a, c a = cos a, sp = sin/3, 
cp = cos/3, tp = tan/?, sp- a = sm(/3 — a), cp_ a = cos(/? — a). 



APPENDIX B: FORM FACTORS 



This appendix lists all the coefficient As in the amplitude of subprocess gg — > ZH, in 
terms of 3- and 4-points one-loop integrals j^j. The diagrams (a) — (e) refer to those in 
Fig. El For convenience, we define abbreviations of one-loop integrals for each diagram as 
following, 

C (a) = C (b) = c{0,s, 0,m 2 f ,m 2 f ,m 2 f ) 

q(c) _ C(m 2 H ,m 2 z , s,m 2 pm 2 ,mj) 

£)( c ) = D(0, m 2 ^, m% } 0, i, s, m 2 , m 2 , m 2 *, m 2 ) 

(j{d) _ C{m 2 Zl m 2 H ,s,m 2 ,m 2 ^m 2 ) 
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D (d) 

C (e) 
D (e) 



D(0, m 2 z , m 2 H , 0, it, s, rrij, m 2 , mj, m 2 ) 

C(0, m 2 H , i, m 2 , m 2 , m 2 ) 

D(m 2 z , 0, m 2 H , 0, it, t, m 2 , m 2 , m 2 , m 2 ) 



Note that in the 
For diagram ( 



4 a) 



following expressions there is an implicit sum over / for / = t, b. 

a) and (b), there are only two coefficients which are not zero, respectively. 



A 



(a) 



2tt 2 



G HZZ( G ffZ _ G fJZ )G f fgiG ff 92 _ 



(a) 



A? 



1 



s — m z 
m f 



2 ^12 



(b) 



5 2tt 2 V l r 1 s-m 2 s + im s T s ° 

(c), the coefficients are the following expressions time an overall factor 



For diagram 
^G^~ H (G{ fZ - G f R fz )G^G^m f 



A, = 4 [<# + Cf + cf - 2D$ + §(D® + £>g>) 



A 2 = 2Cf + 4(Cj CJ + C^) + (t - m l z )Df - 2sD\ 



y(c) 



(c)> 



L>S c) -2^ c) + 4(^ + J D^) 



)(c) 



+ 4(m| + s - i)D[% + 4sL> 



(c) 



)(c) 



13 



^3 
A, 

A, 
A, 
A 7 
As 
A 10 
A 11 
A 12 
A 13 
A u 



-2Ci c) - ACr + s 



(c) 



L>S C) - 2D 2 



4(1© + 



(c) 



-4 



(c) 



,(c) 



(c) 



4 

4 L 

A 9 = -4(^ c) + 2L>g) 
-4(L>{ c) + 2L>g) 
2(D { C) - 2Di c) - AD{f 



2(Z>g + D (c) 



23 



L>i c) + 4L> 



(c) 



2Di c) + 4(Dg + Z> (c) 



D, 



(c) 



2I^ c) 



4(£>£> + £>£' + 



(c) 



23 



-4(L>S c) + 4 C) 



D 



(c) 



Df) - £)W + 2OT + £>g0 



)(c) 



(c) 



The overall factor for diagram (d) is ^G f!H (G T L fZ - G f R fZ )G ff9l G ff92 m f , and the coef- 
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ficients are 



A 2 
A 3 
A, 
A, 
A & 
A 7 

As 
Aw 
A u 
A 12 
A 13 
A u 



A 1 = -2Cf° + 4C ( 2 a> + (m z -u- A)D^> + 2sD\' 



(d) 



+8D$ + 4(u - m 2 z )D ( $ + 4uD^ + 4(u - t)D 



4Cf + (m| 



t) 



D 



(d) 



= 2C [ q ] + AC\ 



(d) 



D, 



(d) 



4 

-4 



o 

d? ■ -ant ■ ». 



2D? 
- 2D? 



12 

j(<0 



(d) 



23 



+ 



4(A2 

l(£>g> + Dg> + £>g') 



D 



(d) 



£>f° + D 



(«0 o/nW 
2 — z l-^13 



4?) 



-4(D$ d) + 2D 

(d) 



(d)^ 

23 , 



= A 9 = -4(£>£" + 2£>g>) 



= 4 



Dl- 



Df 



2(D!J + D 



22 y 



= -2 



D 



(d) 



4(Di d) + D\ a 2 > + D 



(d) 



22 , 



= 2(Dr+4Dg) 



Dj d) - 2D 2 



(d) 



- 4(£>g> + D 2 a 2 > + D 



(d) 



,(d)> 

'23 , 



^(d) _ ^(d) _ ^(d) _ 4(L) 



(d) n (d)N 

12 "I" -^22 / 



4(Df 



2£>g>) 



For diagram (e), the overall factor is ^G ff ~ H ( y G f L fZ -G I R fZ )G f ~ fgi G f ~ f92 m f , and coefficients 



are 



A 1 

A 2 
A 3 
A, 
A 5 
A 6 
A 7 
As 



-2Ct ] + (u- m\)Df + 2(m| + u)Df + AuD ( 2 e) + (m| + s - i^D^ + 8D ( (e) 



oo 



+4m 2 z D$ + 4(m| 



2C,S e) + 8D^ + 2(t 



.(e) 



= 2s 



£)W + 3 £)W + 2(£) w + D w) + 2(2S _ m | )£) w + {i _ m 2 z + 5S)D % 



)(e) 



j(e) 



£>< e) + £#> + 2(£># + £>£') 



Dg) + 4(« - m|)Dg + 4(« - m|) - 5 



23 



(e) 



4sD 



(e) 
23 



(e) 



l(e) 



4(£>< b) + 2D$) 



2Dg; 

(e)> 



-4 
-4 



D( e) + 2(D^+D£) 



(e) 



Dj e) + 3(D\ e> + D 2 e> ) + 2(D\l> + 2Df> + D^) 



>(e) 



l(e) 



(e) 



»(e) 



= 4 



2jW + 2 (Di e 2 ; + £>£') 



(e) 
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A, 
A w 
A ll 
A 12 
A 13 
A 14 



+ Df + 2{D$ + 2D$ + D ( 2 e j 



-2 [L>S e) + 2D! e) + 6I^ e) 



-4 
-2 



i e) +4(D$ + D%+D$+D%) 
~D ( e) + 2( J D 1 e) + Di e) ) + 4( J D 1 e 3 ) + D<$ 
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FIG. 1: Feynman diagrams for the subprocess bb — > ZH 
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FIG. 2: Feynman diagrams for the subprocess gg — > ZH (including only quark loop) 
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FIG. 3: Feynman diagrams for the subprocess gg — ► ZH (including only squark loop) 
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FIG. 4: Total cross sections a (in fb) of pp — > Zh° via bb annihilation (solid lines) compared with 
Drell-Yan ones (dotted lines) at the LHC (a) as functions of m^o for tan/3 = 4, 15 and 40; and (b) 
as functions of tan/3 for m h o = 105 and 115 GeV. 
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FIG. 5: Total cross section a (in fb) of pp — > Zh° via gluon fusion (solid lines) compared with 
complete qq annihilation ones (dotted lines) at the LHC (a) as functions of m h o for tan (3 = 4, 15 
and 40; and (b) as functions of tan/3 for m h o = 105 and 115 GeV. 



16 




100 200 300 400 500 600 700 800 

(a) 



<T[fb] 

10 2 



10 1 - 



10° r- 



10 



-1 _ 



10" 



4n H0 = 400GeV 



10 



(b) 



wz H o[GeV] 




15 20 25 30 35 40 45 50 



tan/3 



FIG. 6: Total cross sections a (in fb) of pp — > via bb annihilation (solid lines) compared 

with Drell-Yan ones (dotted lines) at the LHC (a) as functions of m H o for tan/3 = 4 (curves 
starting from rn ff o=250GeV, 15 (curves starting from m^o=200GeV and 40 (curves starting from 
m#o=100GeV); and (b) as functions of tan/3 for m H o = 200 and 400 GeV. 
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FIG. 7: Total cross section a (in fb) of pp — ► Z-ff via gluon fusion (solid lines) compared with 
complete qq annihilation ones (dotted lines) at the LHC (a) as functions of mjjo for tan (3 = 4, 15 
and 40(the mass ranges are the same as Fig. EJ); and (b) as functions of tan/5 for m#o = 200 and 
400 GeV. 
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FIG. 8: A comparison of the production cross sections of the three Higgs bosons as functions of 
their masses. 
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